Heterostructures coupling ultrathin metal carbides and chalcogenides.

Non-layered transition metal carbides (TMCs) and layered transition metal dichalcogenides (TMDs) are two well-studied material families that have individually received considerable attention over the past century. In recent years, with the shift towards two-dimensional materials and heterostructures, a field has emerged that is focused on the structure and properties of TMC/TMD heterostructures, which through chemical conversion exhibit diverse types of heterostructure configuration that host coupled 2D-3D interfaces, giving rise to exotic properties. In this Review, we highlight experimental and computational efforts to understand the routes to fabricate TMC/TMD heterostructures. Furthermore, we showcase how controlling these heterostructures can lead to emergent electronic transport, optical properties and improved catalytic properties.

Predicted Separation of Acid Gases from Gas Mixtures by Functionalized Porous Aromatic Frameworks.

The selective adsorption of target acid gas molecules from binary gas mixtures by porous aromatic frameworks (PAFs) with two identical functional groups per aromatic ring (PAF-R2) was computationally investigated using grand canonical Monte Carlo simulations. PAF-R2 adsorption was considered for three binary mixtures of small molecular concentrations of acid gas and abundant nitrogen gas (CO2/N2, SO2/N2, and H2S/N2). The results indicate that additional functional groups enhance acid gas loadings and selectivity, compared with pristine PAF and single-functionalized PAFs. Low pressures yield linearly increasing gas loadings and constant selectivity, while high pressures yield much higher adsorption and selectivity. In particular, SO2 loading and selectivity under high pressures are heavily influenced by the PAF's maximum adsorption limit, which can be linked back to the functional groups and their configuration. In summary, PAF-(3,5)-(COOH)2 (nomenclature of PAFs is provided in the Appendix in the Supporting Information) and many other PAF with the same two electron-withdrawing groups are predicted to have great acid gas adsorption and selectivity from gas mixtures, while PAF-(3,5)-(OH)2 (one of PAFs with two identical electron-donating groups) is predicted to have good adsorption and selectivity, especially under elevated pressures. The results of this work can provide insights into various types of PAFs with great selective adsorption ability and their corresponding conditions. The simulation procedures and results may inspire the exploration and screening of other types of PAFs or porous materials, for acid gas absorption.

Single-Step Direct Laser Writing of Multimetal Oxygen Evolution Catalysts from Liquid Precursors.

We investigate a laser direct-write method to synthesize and deposit metastable, mixed transition metal oxides and evaluate their performance as oxygen evolution reaction catalysts. This laser processing method enabled the rapid synthesis of diverse heterogeneous alloy and oxide catalysts directly from cost-effective solution precursors, including catalysts with a high density of nanocrystalline metal alloy inclusions within an amorphous oxide matrix. The nanoscale heterogeneous structures of the synthesized catalysts were consistent with reactive force-field Monte Carlo calculations. By evaluating the impact of varying transition metal oxide composition ratios, we created a stable Fe0.63Co0.19Ni0.18Ox/C catalyst with a Tafel slope of 38.23 mV dec-1 and overpotential of 247 mV, a performance similar to that of IrO2. Synthesized Fe0.63Co0.19Ni0.18Ox/C and Fe0.14Co0.46Ni0.40Ox/C catalysts were experimentally compared in terms of catalytic performance and structural characteristics to determine that higher iron content and a less crystalline structure in the secondary matrix decrease the charge transfer resistance and thus is beneficial for electrocatalytic activity. This conclusion is supported by density-functional theory calculations showing distorted active sites in ternary metal catalysts are key for lowering overpotentials for the oxygen evolution reaction.

Superconductivity enhancement in phase-engineered molybdenum carbide/disulfide vertical heterostructures.

Stacking layers of atomically thin transition-metal carbides and two-dimensional (2D) semiconducting transition-metal dichalcogenides, could lead to nontrivial superconductivity and other unprecedented phenomena yet to be studied. In this work, superconducting α-phase thin molybdenum carbide flakes were first synthesized, and a subsequent sulfurization treatment induced the formation of vertical heterolayer systems consisting of different phases of molybdenum carbide-ranging from α to γ' and γ phases-in conjunction with molybdenum sulfide layers. These transition-metal carbide/disulfide heterostructures exhibited critical superconducting temperatures as high as 6 K, higher than that of the starting single-phased α-Mo2C (4 K). We analyzed possible interface configurations to explain the observed moiré patterns resulting from the vertical heterostacks. Our density-functional theory (DFT) calculations indicate that epitaxial strain and moiré patterns lead to a higher interfacial density of states, which favors superconductivity. Such engineered heterostructures might allow the coupling of superconductivity to the topologically nontrivial surface states featured by transition-metal carbide phases composing these heterostructures potentially leading to unconventional superconductivity. Moreover, we envisage that our approach could also be generalized to other metal carbide and nitride systems that could exhibit high-temperature superconductivity.

Optimized utilization of COMB3 reactive potentials in LAMMPS.

An investigation to optimize the application of the third-generation charge optimized many-body (COMB3) interatomic potential and associated input parameters was carried out through the study of solid-liquid interactions in classical molecular dynamics simulations. The rates of these molecular interactions are understood through the wetting rates of water nano-droplets on a bare copper (111) surface. Implementing the Langevin thermostat, the influence of simulation time step, the number of atoms in the system, the frequency at which charge equilibration is performed, and the temperature relaxation rate are all examined. The results indicate that time steps of 0.4 fs are possible when using longer relaxation times for the system temperature, which is almost double the typical time step used for reactive potentials. The use of the charge equilibration allows for a fewer atomic layers to be used in the Cu slab. In addition, charge equilibrium schemes do not need to be performed every time step to ensure accurate charge transfer. Interestingly, the rate of wetting for the nanodroplets is dominantly dependent on the temperature relaxation time, which is predicted to significantly change the viscosity of the water droplets. This work provides a pathway for optimizing simulations using the COMB3 reactive interatomic potential.

Effects of surface charge and cluster size on the electrochemical dissolution of platinum nanoparticles using COMB3 and continuum electrolyte models.

We study the site-dependent dissolution of platinum nanoparticles under electrochemical conditions to assess their thermodynamic stability as a function of shape and size using empirical molecular dynamics and electronic-structure models. The third-generation charge optimized many-body potential is employed to determine the validity of uniform spherical representations of the nanoparticles in predicting dissolution potentials (the Kelvin model). To understand the early stages of catalyst dissolution, implicit solvation techniques based on the self-consistent continuum solvation method are applied. It is demonstrated that interfacial charge and polarization can shift the dissolution energies by amounts on the order of 0.74 eV depending on the surface site and nanoparticle shape, leading to the unexpected preferential removal of platinum cations from highly coordinated sites in some cases.

The structure of graphene on graphene/C60/Cu interfaces: a molecular dynamics study.

Two experimental studies reported the spontaneous formation of amorphous and crystalline structures of C60 molecules intercalated between graphene and a surface. The findings observed included interesting phenomena ranging from reaction between fullerene C60s ('C60s' stands for plural of C60) under graphene to graphene sheets sagging between C60s and control of strain in these sheets. Motivated by this work, we performed fully atomistic reactive molecular dynamics simulations to investigate the formation and thermal stability of graphene sheet wrinkles as well as graphene attachment to and detachment from a surface when the sheet is laid over a previously distributed array of C60 molecules on a copper surface at different temperatures. As graphene compresses the C60s against the surface, and graphene attachment to the surface in between C60s depends on the height of the wrinkles in the graphene sheet, configurations with both frozen and non-frozen fullerenes were investigated in the simulations in order to examine the experimental result of stable, sagged graphene sheets when the distance between C60s is about 4 nm and the height of the wrinkles in the sheet is about 0.8 nm. Below a distance of 4 nm between fullerenes, the graphene is predicted to become locally suspended and less strained. The simulations predict that this happens when the fullerenes can deform under the compressive action of the graphene sheet. If the fullerenes are kept frozen, spontaneous 'blanketing' of graphene is predicted only when the distance between neighbouring C60s is equal to or great than about 7 nm. These predictions agree with a mechanical model relating the rigidity of a graphene sheet to the energy of graphene-surface adhesion. This work further reveals the structure of intercalated molecules and the role of stability and sheet wrinkling on the preferred configuration of graphene. This study thus might assist in the development of two-dimensional confined nanoreactors for chemical reactions.

Multi-Step Topochemical Pathway to Metastable Mo2AlB2 and Related Two-Dimensional Nanosheet Heterostructures.

The rational synthesis of metastable inorganic solids, which is a grand challenge in solid-state chemistry, requires the development of kinetically controlled reaction pathways. Topotactic strategies can achieve this goal by chemically modifying reactive components of a parent structure under mild conditions to produce a closely related analogue that has otherwise inaccessible structures and/or compositions. Refractory materials, such as transition metal borides, are difficult to structurally manipulate at low temperatures because they generally are chemically inert and held together by strong covalent bonds. Here, we report a multistep low-temperature topotactic pathway to bulk-scale Mo2AlB2, which is a metastable phase that has been predicted to be the precursor needed to access a synthetically elusive family of 2-D metal boride (MBene) nanosheets. Room-temperature chemical deintercalation of Al from the stable compound MoAlB (synthesized as a bulk powder at 1400 °C) formed highly strained and destabilized MoAl1-xB, which was size-selectively precipitated to isolate the most reactive submicron grains and then annealed at 600 °C to deintercalate additional Al and crystallize Mo2AlB2. Further heating resulted in topotactic decomposition into bulk-scale Mo2AlB2-AlOx nanolaminates that contain Mo2AlB2 nanosheets with thickness of 1-3 nm interleaved by 1-3 nm of amorphous aluminum oxide. The combination of chemical destabilization, size-selective precipitation, and low-temperature annealing provides a potentially generalizable kinetic pathway to metastable variants of refractory compounds, including bulk Mo2AlB2 and Mo2AlB2-AlOx nanosheet heterostructures, and opens the door to other previously elusive 2-D materials such as 2-D MoB (MBene).

Carbon doping of WS2 monolayers: Bandgap reduction and p-type doping transport.

Chemical doping constitutes an effective route to alter the electronic, chemical, and optical properties of two-dimensional transition metal dichalcogenides (2D-TMDs). We used a plasma-assisted method to introduce carbon-hydrogen (CH) units into WS2 monolayers. We found CH-groups to be the most stable dopant to introduce carbon into WS2, which led to a reduction of the optical bandgap from 1.98 to 1.83 eV, as revealed by photoluminescence spectroscopy. Aberration corrected high-resolution scanning transmission electron microscopy (AC-HRSTEM) observations in conjunction with first-principle calculations confirm that CH-groups incorporate into S vacancies within WS2. According to our electronic transport measurements, undoped WS2 exhibits a unipolar n-type conduction. Nevertheless, the CH-WS2 monolayers show the emergence of a p-branch and gradually become entirely p-type, as the carbon doping level increases. Therefore, CH-groups embedded into the WS2 lattice tailor its electronic and optical characteristics. This route could be used to dope other 2D-TMDs for more efficient electronic devices.

Three-dimensional atomic scale electron density reconstruction of octahedral tilt epitaxy in functional perovskites.

Octahedral tilts are the most ubiquitous distortions in perovskite-related structures that can dramatically influence ferroelectric, magnetic, and electronic properties; yet the paradigm of tilt epitaxy in thin films is barely explored. Non-destructively characterizing such epitaxy in three-dimensions for low symmetry complex tilt systems composed of light anions is a formidable challenge. Here we demonstrate that the interfacial tilt epitaxy can transform ultrathin calcium titanate, a non-polar earth-abundant mineral, into high-temperature polar oxides that last above 900 K. The comprehensive picture of octahedral tilts and polar distortions is revealed by reconstructing the three-dimensional electron density maps across film-substrate interfaces with atomic resolution using coherent Bragg rod analysis. The results are complemented with aberration-corrected transmission electron microscopy, film superstructure reflections, and are in excellent agreement with density functional theory. The study could serve as a broader template for non-destructive, three-dimensional atomic resolution probing of complex low symmetry functional interfaces.

Nanoscale Structure and Dynamics of Water on Pt and Cu Surfaces from MD Simulations.

The interaction of liquid water with Pt(111) is investigated with classical molecular dynamics (MD) simulations, where the forces are determined using the third-generation charge optimized many-body (COMB3) interatomic potential. In cases of sub-monolayer water coverage, the parameterized empirical potential predicts experimentally observed and energetically favorable √37 and √39 reconstructed water structures with "575757" di-interstitial defects. At both sub-monolayer and multilayer water coverages, the structure of the first wetting layer of liquid water on Pt(111) exhibits a characteristic distribution where the molecules form two distinct buckled layers as a result of the interplay between water-metal adsorption and water-water hydrogen bonds. The dynamic spreading rate of water nanodroplets on large Pt surfaces (>200 nm2) characterized by molecular kinetic spreading theory is an order of magnitude slower than the molecular kinetic rate of the same droplet on close-packed Cu surfaces due to variation in molecular distributions at the water-metal interface. These nanoscale MD simulation predictions using the COMB3 interatomic potential demonstrate the capability of capturing both many-body interactions between H2O and Pt or Cu and hydrogen bonding in liquid water.

Applied Potentials in Variable-Charge Reactive Force Fields for Electrochemical Systems.

An atomic description of water dynamics and electrochemical properties at electrode-electrolyte interfaces is presented using molecular dynamics with the third generation of the charge-optimized many-body (COMB3) potential framework. Externally applied potentials in electrochemical applications were simulated by offsetting electronegativity on electrode atoms. This approach is incorporated into the variable charge scheme within COMB3 and is used to investigate electrochemical systems consisting of two Cu electrodes and a water electrolyte with varying concentrations of hydroxyls (OH-) and protons (H+). The interactions between the electronegativity offset method and the charge equilibration method in a variable charge scheme are analyzed. In addition, a charge equilibration method for electrochemical applications is proposed, where the externally applied potentials are treated by the electronegativity offset on the electrodes thus enforcing charge neutrality on the electrolyte. This method is able to qualitatively capture the relevant electrochemistry and predict consistently correct voltages with precalibration.

Computational Study of Low Interlayer Friction in Tin+1Cn (n = 1, 2, and 3) MXene.

The friction of adjacent Tin+1Cn (n = 1, 2, and 3) MXene layers is investigated using density functional theory (DFT) calculations and classical molecular dynamics simulations with ReaxFF potentials. The calculations reveal the sliding pathways in all three MXene systems with low energy barriers. The friction coefficients for interlayer sliding are evaluated using static calculations. Both DFT and ReaxFF methods predict friction coefficients between 0.24 and 0.27 for normal loads less than 1.2 GPa. The effect of titanium (Ti) vacancies in sublayers and terminal oxygen (O) vacancies at surfaces on the interlayer friction is further investigated using the ReaxFF potential. These defects are found to increase the friction coefficients by increasing surface roughness and creating additional attractive forces between adjacent layers. However, these defective MXenes still maintain friction coefficients below 0.31. We also consider functionalized Ti3C2 MXene terminated with -OH and -OCH3 and find that compared to the -O-terminated surface both groups further reduce the interlayer friction coefficient to 0.10-0.14.

Graphene-Titanium Interfaces from Molecular Dynamics Simulations.

Unraveling the physical and chemical properties of graphene-metal contacts is a key step toward the development of graphitic electronic nanodevices. Although many studies have revealed the way that various metals interact with graphene, few have described the structure and behavior of large pieces of graphene-metal nanostructures under different conditions. Here, we present the first classical molecular dynamics study of graphene-titanium (G-Ti) structures, with and without substrates. Physical and chemical properties of equilibrium structures of G-Ti interfaces with different amounts of titanium coverage are investigated. Adhesion of Ti films on graphene is shown to be enhanced by the vacancies in graphene or the electrostatic influence of substrates. The dynamics of pristine G-Ti structures at different temperatures on planar and nonplanar substrates are investigated, and the results show that G-Ti interfaces are thermally stable, that is, not prone to any reaction toward the formation of titanium carbide.

Two-Dimensional Intrinsic Half-Metals With Large Spin Gaps.

Through a systematic search of all layered bulk compounds combined with density functional calculations employing hybrid exchange-correlation functionals, we predict a family of three magnetic two-dimensional (2D) materials with half-metallic band structures. The 2D materials, FeCl2, FeBr2, and FeI2, are all sufficiently stable to be exfoliated from bulk layered compounds. The Fe2+ ions in these materials are in a high-spin octahedral d6 configuration leading to a large magnetic moment of 4 µB. Calculations of the magnetic anisotropy show an easy-plane for the magnetic moment. A classical XY model with nearest neighbor coupling estimates critical temperatures, Tc, for the Berezinskii-Kosterlitz-Thouless transition ranging from 122 K for FeI2 to 210 K for FeBr2. The quantum confinement of these 2D materials results in unusually large spin gaps, ranging from 4.0 eV for FeI2 to 6.4 eV for FeCl2, which should defend against spin current leakage even at small device length scales. Their purely spin-polarized currents and dispersive interlayer interactions should make these materials useful for 2D spin valves and other spintronic applications.

Topology-Scaling Identification of Layered Solids and Stable Exfoliated 2D Materials.

The Materials Project crystal structure database has been searched for materials possessing layered motifs in their crystal structures using a topology-scaling algorithm. The algorithm identifies and measures the sizes of bonded atomic clusters in a structure's unit cell, and determines their scaling with cell size. The search yielded 826 stable layered materials that are considered as candidates for the formation of two-dimensional monolayers via exfoliation. Density-functional theory was used to calculate the exfoliation energy of each material and 680 monolayers emerge with exfoliation energies below those of already-existent two-dimensional materials. The crystal structures of these two-dimensional materials provide templates for future theoretical searches of stable two-dimensional materials. The optimized structures and other calculated data for all 826 monolayers are provided at our database (https://materialsweb.org).

Effect of Surface Chemistry on Water Interaction with Cu(111).

The interfacial dynamics of water in contact with bare, oxidized, and hydroxylated copper surfaces are examined using classical molecular dynamics (MD) simulations. A third-generation charge-optimized many-body (COMB3) potential is used in the MD simulations to investigate the adsorption of water molecules on Cu(111), and the results are compared to the findings of density functional theory (DFT) calculations. The adsorption energies and structures predicted by COMB3 are generally consistent with those determined with DFT. The COMB3 potential is then used to investigate the wetting behavior of water nanodroplets on Cu(111) at 20, 130, and 300 K. At room temperature, the simulations predict that the spreading rate of the base radius, R0, of a water droplet with a diameter of about 1.5 nm exhibits a spreading rate of R0 ≈ t(0.16) and a final base radius of 3.5 nm. At 20 and 130 K, water droplets are predicted to retain their structure after adsorption on Cu(111) and to undergo minimal spreading in agreement with scanning tunneling microscopy data. When the same water droplet encounters a reconstructed, oxidized Cu(111) surface, the classical MD simulations predict wetting with a spreading rate of R ≈ t(0.14) and a final base radius of 3.0 nm. Similarly, our MD simulations predict a spreading rate of R ≈ t(0.14) and a final base radius of 2.5 nm when water encounters OH-covered Cu(111). These results indicate that oxidation and hydroxylation cause a reduction in the degree of spreading and final base radius that is directly associated with a decreased spreading rate for water nanodroplets on copper.

Mapping Chemical Selection Pathways for Designing Multicomponent Alloys: an informatics framework for materials design.

A data driven methodology is developed for tracking the collective influence of the multiple attributes of alloying elements on both thermodynamic and mechanical properties of metal alloys. Cobalt-based superalloys are used as a template to demonstrate the approach. By mapping the high dimensional nature of the systematics of elemental data embedded in the periodic table into the form of a network graph, one can guide targeted first principles calculations that identify the influence of specific elements on phase stability, crystal structure and elastic properties. This provides a fundamentally new means to rapidly identify new stable alloy chemistries with enhanced high temperature properties. The resulting visualization scheme exhibits the grouping and proximity of elements based on their impact on the properties of intermetallic alloys. Unlike the periodic table however, the distance between neighboring elements uncovers relationships in a complex high dimensional information space that would not have been easily seen otherwise. The predictions of the methodology are found to be consistent with reported experimental and theoretical studies. The informatics based methodology presented in this study can be generalized to a framework for data analysis and knowledge discovery that can be applied to many material systems and recreated for different design objectives.

Charge optimized many-body (COMB) potential for dynamical simulation of Ni-Al phases.

An interatomic potential for the Ni-Al system is presented within the third-generation charge optimized many-body (COMB3) formalism. The potential has been optimized for Ni3Al, or the γ' phase in Ni-based superalloys. The formation energies predicted for other Ni-Al phases are in reasonable agreement with first-principles results. The potential further predicts good mechanical properties for Ni3Al, which includes the values of the complex stacking fault (CSF) and the anti-phase boundary (APB) energies for the (1 1 1) and (1 0 0) planes. It is also used to investigate dislocation propagation across the Ni3Al (1 1 0)-Ni (1 1 0) interface, and the results are consistent with simulation results reported in the literature. The potential is further used in combination with a recent COMB3 potential for Al2O3 to investigate the Ni3Al (1 1 1)-Al2O3 (0 0 01) interface, which has not been modeled previously at the classical atomistic level due to the lack of a reactive potential to describe both Ni3Al and Al2O3 as well as interactions between them. The calculated work of adhesion for this interface is predicted to be 1.85 J m(-2), which is in agreement with available experimental data. The predicted interlayer distance is further consistent with the available first-principles results for Ni (1 1 1)-Al2O3 (0 0 0 1).